package scu.maqiang.matlabcases;

import scu.maqiang.fes.BVPType;
import scu.maqiang.mesh.Mesh2T3;
import scu.maqiang.mesh.MeshFormat;
import scu.maqiang.possion.Possion2DT31;

public class PoissonCase1 {

	public static void main(String[] args) {
		String fileName = "matlabMeshes\\PoissonCase1Mesh.msh";
		Mesh2T3 mesh2D = new Mesh2T3().readMesh(fileName, MeshFormat.msh);
		Possion2DT31 p2DT31 = new Possion2DT31(mesh2D);
		int[] DirichletBoundaryLabel = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};
		int[] NeumannBoundaryLabel = {13, 14, 15, 16, 17, 18, 19};
		p2DT31.GeneralPossionProblem(BVPType.COMMON, 1.0, 10,
				                     DirichletBoundaryLabel, 0.0, 
				                     NeumannBoundaryLabel, -5, 
				                     null, 0.0, 0.0, "matlabResults\\PoissonCase1Result.dat");
		
	}

}
